Summary 


AeroPropulsoServoElasticity: Dynamic Modeling of 
the Variable Cycle Propulsion System 


This presentation was made at the 2012 Fundamental Aeronautics Program Technical 
Conference and it covers research work for the Dynamic Modeling of the Variable 
cycle Propulsion System that was done under the Supersonics Project, in the area of 
AeroPropulsoServoElasticity. 

The presentation covers the objective for the propulsion system dynamic modeling 
work, followed by the work that has been done so far to model the variable Cycle 
Engine, modeling of the inlet, the nozzle, the modeling that has been done to model 
the affects of flow distortion, and finally presenting some concluding remarks and 
future plans. 
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Outline 


* Objective 

* Variable Cycle Engine (VCE) Propulsion System Modeling 

- Engine Modeling 

- Inlet modeling 

- Nozzle Modeling 

* Parallel Flow Path Modeling 



Concluding Remarks/Future 
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Objective 

What is AeroPropulsoServoElasticity (APSE), Why? 




Vehicle ASE Model 

Propulsion System Thrust 




> Integrated 
Modeling 
&Controls Design 


> Vehicle Stability 


> Ride quality 



> Cruise Efficiency 


Integrated APSE Model 
(NASA LaRC in collaboration with NASA GRC) 


Design and Analysis 



Key Milestone Progress 

N+2/N+3 P ropulsion System Modeling 



Objective 

• Develop concept Variable Cycle Engine (VCE) 
dynamic propulsion system model 

• Develop propulsion system operation control designs 

Approach 

• Modify prior developed engine models & develop 
additional component models 

• Develop CFD inlet and nozzle models w/ variable 
geometries 

• Generate linear models and feedback control designs 

• Generate engine operation control schedules 

Significance 

• Provides a valuable platform for controls design to 
improve performance & disturbance suppression 

• Provides a platform for dynamically verifying 
propulsion system concept designs 

• Provides propulsion system for ASE integration and 
for overall controls & vehicle performance studies 
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VCE Component Model Diagram 


Variable Cycle Engine Model 


• Dual Spool variable cycle - High bypass at low altitudes to low bypass high 
altitudes 


• Noise abatement for overland flight 
-- Through external bypass & through nozzle design 


• Cycle analysis conducted in NPSS - provided geometries and component 
performance characteristics for dynamic model 


main bypass external bypass 



core flow 
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Component Modeling - Roadmap & Approach 


Development Roadmap 

Original component models developed 
based on J85-13 engine 

Many of J85-13 component models and 
methods directly utilized for VCE w/ the 
appropriate maps and geometries 

Some new component models 
developed (ducts, mixers, splinters, 
dual core) -VCEV.l 
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Volume Dynamics 


For some components need to develop 
detailed models - like CFD for inlet & 
nozzles 


Need to develop fully operational 
engine (control schedules) - 
Methodology developed w/ J85-13 

Parallel flow paths for distortion & 
boundary layer effects 

Propulsion & ASE integration - 
Interfaces, controls and performance 
studies 


Continuity of mass, momentum & energy 
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VCESplitter 


Variable Cycle Engine Model Components 
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Bleed HPT2 
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Bleed HPT1 
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VCE Engine Results 
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Initial objective is VCE model development 

• Control design effort light; hold model together 
-- But designed for higher bandwidth controls 

for disturbance attenuation 

• Engine has higher response capability of ~ 70 
rad/sec on high side (~40 rad/sec typically used) 

• Potential to use higher response capability to 
design for better disturbance attenuation, safety 
margins, and engine efficiency 
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VCE Engine Speed and Thrust 
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time, (s) 




3 3.1 3.2 3.3 3.4 

time, (s) 


• Nominal VCE propulsion system thrust 
44,100 Nor 9,914 Ibf 

• A 1% change in fan speed causes 
2.9% change in thrust 

• Thrust response more underdamped - design 
of speed controller also needs to consider 
thrust response 


VCE Engine Atmospheric Disturbance and Thrust 


Thrust response w/ Atmospheric Disturbance 

With no external compression inlet & no ID 
CFD for nozzles 

• Case 1; eddy dissipation rate 4x average of 
North Atlantic cruise altitudes; integral length 
scale typical (equivalent to atmospheric 
turbulence patch size of ~ 11 km); max locally 
dissipating wind speeds 80 mph 

-- Results in thrust variations up to ~ 

5000 N or 1124 lb 

• Case 2; eddy dissipation rate worst 
recorded; integral length scale typical; max 
dissipating wind speeds 150 mph 

-- Results in thrust variation up to ~ 

9000N or 20241b 
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Variable Cycle Propulsion System Studies 

Preliminary - Thrust Spectral for Coupling to AeroServoElastic (ASE) Modes 


• Study based on V.l initial 
variable cycle engine 
modeling 


• Atmospheric turbulence 
model w/ eddy dissipation 
rates & momentary wind 
gusts up to 180 mph 


• Study shows potentially 
significant trust dynamics 
to warrant detailed APSE 
modeling and analysis 
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Supersonic Inlets Modeling 



Started with Mixed Compression Supersonic inlets 
(results presented last year) 


'External 
Supersonic 
Compression 


Throat / terminal 
shock 


Internal 

Supersonic 

Diffuser 


Subsonic 

Diffuser ^ n 8' ne 
Face 


Now focusing on external compression axisymmetric Inlets 
-- Better overall performance for Mach 1.8 or less 



External Compressi on Inlet Modeling - Approach 

Computational Domain 

A. l-D compressible flow cells w/ flow propagation delay dynamics and averaging 
flows at shock boundary 

B. Quasi l-D CFD compressible flow cells w/ leakage fluxes estimation and artificial 
viscosity 

C. Quasi l-D CFD compressible flow cells w/ artificial viscosity 



A-B. Moving computational domains 
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External Compression Inlet Modeling - Challenges 


Challenges 

- Developing generalized formulations for conservation flux 
leakages across sonic boundary - Method hasn't worked yet 

- Sensing the shock position to switch between compressible 
flow cells and quasi ID CFD cells - Moving Domain 

- Determined mass flow leakage based on test data for various 
engine face back pressures to calculate leakage fluxes - 
Approach worked but is not generalized 

- Remaining issue for inlet dynamics Conical 
compressible flow field inherently 2D and 
3D for pitch variations 
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External Compression Inlet Results 
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Comparing test and Simulation Results 
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Pressure vs. Inlet Domain 



Pressure profile by ramping back pressure 


Difference In Shock Position 


Back Pressure 
(N/m 2 ) 

Test Data Shock Position 
(Cell) 

Simulation Shock Position 
(Cell) 
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Nozzle Modeling 



Obiective/Approach 

• Develop ID CFD model for exit nozzles for thrust dynamics (before used nozzle 
lump volume and chocked compressible flow function) 

-- Chosen method: MacCormack's predictor-corrector technique assuming 
subsonic-supersonic isentropic nozzle flow 


Step one - develop model for generic Convergent-Divergent (CD) nozzle 
geometry - Simple shape profile w/ actual N+3 throat and areas 

Step two - develop model for more complex supersonic engine- 
nozzle concept geometry 


Inlet 
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CFD Method- Predictor Step 
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Static temperature, (K) normalized area & density 


Results 

(so far steady state - no freq responses) 




nozzle length (m) 


• Generic model verified against 
results reported in literature 



nozzle length (m) 
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Parallel Compressor Modeling 


Objective 

Develop parallel component flow path models to study effect of 
distortion on propulsion system dynamics and APSE 

• New model derived in 

cylindrical coordinates - Euler 


• Allows modeling of disturbance 

from changing flight conditions Original model 

(pitch, yaw, roll, etc) Stage-by-stage, single flow path 



• Inlet conditions of Pressure, 
Temperature & outlet 
conditions of mass flow rate 

• Path ratio of /?i - scaling 
mass flow rate of stage 
maps by path ratio 
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New Model - Multiple Interacting Flow Paths 


Parallel Compressor Modeling Approach 

Conservation Dynamics in 2D Cylindrical Coordinates 

□ Equations were derived in cylindrical coordinates for compressible & 
inviscid flow, assuming flow properties do not vary in the radial direction 



Conservation Equations 
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Pressure ratio 


Parallel Compressor Modeling Results 


Pressure distortion of approximately 0.1% piji 
applied to path 1 


P4.. 

P3.. 


P2, T2 






mass flow rate 

Pressure disturbance moves Path 1, Path 3 operating points to surge line 


• Would experience cascading stall if mass flow rate was not held constant 
(as with engine) 


Pressure disturbance (pa) 


Parallel Compressor Modeling Results 



Square wave distortion applied to compressor input, path 1 
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Parallel Compressor Modeling Results 


1 st Stage at 100% Speed w/ 1300pa (0.16%) Distortion on Sector 2 & 4 


Stage Rotational Velocity Stallva.12 



Stall Pattern - From Back to Front of Compressor 

(0 Normal, >0 Stall) 
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Conclusion/Future Plans 

• Developed first version of VCE model and preliminary analysis 

• Develop complete dynamic VCE propulsion system models and 
control designs (feedback and operation schedules) 

• Develop Integrated APSE system models, integrated vehicle 
controls, and conduct APSE studies 

• Close integration between NPSS and APSE (already started) 

Additional Possibilities of this Research 

• Integrate w/ NPSS to develop a complete cycle deck design and 
verification package and controls development platform/Rig 

• With gas dynamic model explore higher bandwidth controls to 
reduce stall margins and improve efficiency and design advanced 
controls to improve flight safety and operability 
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